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ABSTRACT 

We have used the Very Long Baseline Array to image 18 quasars with foreground damped 
Lyman-a systems (DLAs) at 327, 610 or 1420 MHz, to measure the covering factor / of 
each DLA at or near its redshifted Hi 21cm line frequency. Including six systems from the 
literature, we find that none of 24 DLAs at 0.09 < z < 3.45 has an exceptionally low covering 
factor, with / — 0.45 - 1 for the 14 DLAs at z > 1 .5, / ~ 0.41 - 1 for the 10 systems at z < 1, 
and consistent covering factor distributions in the two sub-samples. The observed paucity of 
detections of Hi 21cm absorption in high-z DLAs thus cannot be explained by low covering 
factors and is instead likely to arise due to a larger fraction of warm Hi in these absorbers. 

Key words: quasars: individual : quasars: images- galaxies: ISM 



1 INTRODUCTION 

Damped Lyman-a systems (DLAs), selected on the basis of their 
high Hi column densities (A'hi > 2 x 10 2l) cm~ 2 ) in quasar absorp- 
tion spectra, have long b een identified as the progenitors of nor- 
mal present-day galaxies dWolfe et al. 1986), and have hence been 
the subject of much research. Despite this, the nature of high-z 
DLAs, physical conditions in them, and th eir evolution with red- 
shift, are still matters of dispute today (e.g. lWolfe et alj|2005h . An 
issue of recent controversy is the temperature distribution of Hi 
in DLAs. For radio-loud background quasars, a comparison be- 
tween the Hi column density and the optical depth in the redshifted 
Hi 21cm line yields the DLA spin temperature (T s ). For optically- 
thin multi-phase absorption, T s is the column-density-weighted har- 
monic mean of the spin temperatures of the different phases along 
the line of sight. For nearly three decades, the spin temperatures of 
high-z DLAs have been found to be systematically hig her than val- 
ues seen in the Milky Way or local spiral disks (e.g.lWolfe & Davisl 
1 19791 : ICarilli et all Il99rj iKanekar & Chengaludl2003l) ; the first 
low-7^ DLA at high redshift was only detected very recently 
dYork et al] [2007). The simplest interpretation of the observations 
is that Hi in high-z DLAs is predominantly warm, unlike the situa- 
tion in the Galaxy dCarilli et al.ll996r . lChengarur & Kanekai|2000l : 
iKanekar & Chengaluj|2003h . Conversely. IWolfe et all J20031) used 
CII* absorption lines to argue that roughly half the DLAs at z > 2 
contain significant fractions of cold Hi. 

A problem in measuring DLA spin temperatures is that low- 
frequency radio emission is often extended over large angular 
scales, far larger than the size of a galaxy. Redshifted Hi 21cm ab- 
sorption studies are usually carried out with single dishes or short- 
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baseline interferometers, of fairly poor angular resolution, > 10", 
or > 80 kpc at z ~ 3. It is hence often not certain that the quasar 
radio emission is entirely covered by the foreground DLA. If some 
fraction of the quasar emission leaks out around the DLA, the in- 
ferred spin temperature will be an over-estimate. The fraction of the 
radio emissi on occulted by the DL A is referred to as the covering 
factor / (e.g. lBriggs & Wolfell 19831) . 

Recently, ICurran et al. I d2005h have emphasized the problem 
of unknown covering facto rs in high-z DLAs, pointing out that the 
high estimated T s values o flKanekar & Chen galur ( 2003) (hereafter 
KC03) could merely stem from very low DLA covering factors, 
f « l. lCurran & Webbl ( feOOfj) also argued that low DLA covering 
factors at high redshifts arise due to a geometrical effect, as the sim- 
ilarity in angular diameter distances of the high-z DLAs and their 
background quasars implies that a high-z DLA is "les s effective" at 
coveri ng the quasar than a low-z absorber. Note that ICurran et al.l 
(2005) find that quasars of angular extent < 0.1" are adequately 
covered even by compact foreground DLAs. 

The most direct way of resolving this issue is by high spatial 
resolution very long baseline interferometry (VLBI) studies in th e 
redshifted Hi 21cm line (e.g. lBriggs et alj| 19891 ; lEane et alj|2000l) . 
Unfortunately, the low-frequency coverage of current VLBI facili- 
ties is quite limited, implying that such observations are only pos- 
sible for very few sources. Alternatively, one could measure the 
compact flux density arising from the quasar core in VLBI contin- 
uum images at or near the redshifted Hi 21cm line frequency and 
then estimate the DLA covering factor by comparing the VLBI flux 
density with the flux density measured with short baseline interfer- 
omete rs or single dishes (e.g. Brigg s &Woifelll983l; l Kanek ar et al.l 
l2007h . This is the approach that we follow here, based on Very 
Long Baseline Array (VLB A) imaging of 18 quasars with fore- 
ground DLAs at z ~ 0.24 - 3.45. 
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Figure 1. VLBA images of the 18 quasars of the sample, ordered (from top left) by DLA redshift. The map frequency is at the top right of each panel. 
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2 OBSERVATIONS, DATA ANALYSIS AND RESULTS 

The VLBA 1.4 GHz, 610 MHz and 327 MHz receivers were used 
to observe the 18 quasars of our sample between 2002 March and 
2006 June (projects BK89 and BK131), with observing frequencies 
close to the redshifted Hi 21cm line frequency of the foreground 
DLA. Three quasars with DLAs at z ~ 0.24 — 0.44 were observed 
at 1.4 GHz, two, with DLAs at z ~ 0.52 - 0.86, at 610 MHz, and 
13, with DLAs at z ~ 1.78 - 3.45, at 327 MHz. Total bandwidths 
of 4, 8 and 16 MHz were used for the observations at 610 MHz, 
327 MHz and 1420 MHz, respectively, with 2-bit sampling, two 
polarizations, 16 spectral channels and on-source integration times 
of 1 - 4 hours. Strong fringe finders (3C345, 3C454.3 or 3C286) 
were observed every few hours, to calibrate the shape of the pass- 
band. Phase-referencing was not used. The U-V coverage was poor 
in some cases, as all the VLBA antennas were not available, result- 
ing in asymmetric or larger synthesized beams. 

The data were analysed in "classic" AIPS, using standard 
techniques. Initial procedures included ionospheric corrections, 
editing of radio frequency interference, amplitude calibration (us- 
ing the measured antenna gains and system temperatures), pass- 
band calibration and fringe-fitting for the delay rates, before the 
data for each target were averaged in frequency to a single-channel 
dataset. For each source, this was followed by a number of cycles 
of self-calibration and imaging to determine the antenna gains, un- 
til no improvement was seen on further self-calibration. The final 
images obtained from the above procedure are shown in Fig. [T] in 
order of increasing redshift, with the synthesized beams listed in 
column (6) of Table [T] The root-mean-square noise in off-source 
regions in the images was measured to be ~ 0.5 — 3 mJy/Bm. 

The task UVFIT was used to fit elliptical gaussian models 
to the visibility data, to measure the compact flux density. Most 
sources in Fig.Q]are core-dominated, with only a small fraction of 
the flux density in extended structure. For some sources with visible 
extensions (e.g. 1229-021, 1243-072 and 0201+113), a single- 
component model was found sufficient to recover the "cleaned" 
flux density. The only exceptions are 0952+179 (z abs ~ 0.238), 
1157+014 (z abs ~ 1.944) and 0438-436 (z abs ~ 2.347), where 
all the "cleaned" flux density could be recovered only with a 2- 
component model. 

Table Q] summarizes the results of the VLBA observations, 
grouping the sources by observing band, in order of increasing 
redshift. For each source, the first four columns contain the ob- 
serving band, source name, DLA redshift and redshifted Hi 21cm 
frequency. Column (5) lists the total flux densities at the VLBA ob- 
serving frequency. For the three 1 .4 GHz targ ets, this was obtaine d 
from the 1.4 GHz NRAO VLA Sky Survey dCondon et alii 19981) . 
For most 327 MHz targets (11/13), we use the flux density mea- 
sured at nearby freque ncies [e.g. at 327 MHz , from the Westerbork 
Northern Sky Survey jRengelink et alll 19971), 365 MH z, from the 
Texas survey jDouglas et al.ll 19961) , ICarilli et al.l i 19961) or KC03]. 
Four sources (0827+243, 0454+039, 0405-331 and 0438-436) do 
not have measurements in the literature near the VLBA observing 
frequencies; their flux densities were hence estimated from their 
low-frequency (< 1 GHz) spectral indices. 

Columns (6 — 10) contain (6) the VLBA synthesized beam, 
(7) the flux density obtained on fitting a gaussian model with 
UVFIT, (8) the deconvolved angular sizes of the gaussian com- 
ponents, (9) the corresponding spatial extents of these compo- 
nents, at the redshift of the foreground DLA [we use (£l\,Q. m ,h) = 
(0.7,0.3,0.7) in this paper] and (10) the covering factor of the com- 
pact radio emission at the observing frequency, obtained by divid- 



ing the "core" flux density by the total flux density. In three cases 
of two source components, the one that is more compact is identi- 
fied with the core. We emphasize that the deconvolved sizes listed 
in columns (8) and (9) are upper limits, due to the possibility of 
residual phase errors in the data (e.g. due to fluctuations on time 
scales shorter than the self-calibration interval). This is especially 
true for the 327 MHz and 610 MHz results. 



3 OTHER SOURCES 

A few other quasars with foreground DLAs have estimates of the 
covering factor from VLBI studies. These are discuss e d belo w. 

(1) 0738+313 (Zab s ~ 0.0912, 0.2212): lLane et all bOOd) used 
1302 MHz VLBA observations to estimate that the covering fac- 
tor of the foreground DLAs is / ~ 0.98. 

(2) 1127-145 (z abs ~ 0.3127): lBondi et al J fl996l) measured a flux 
density of ~ 5.0 Jy in their 1.6 GHz VLBI image, compared to a 
single-dish flux density of ~ 5 .6 Jy. This gives / ~ 0.89. 

(3) 0235+164 (z abs ~ 0.524): IWolfe et al J dl978l) used 931 MHz 
VLBI observations to s how that the siz e of the quasar core is 
< 6 mas. More recently. iFrev et al J 120001) found that the 1.6 GHz 
VLBI core flux density (with a sub-mas beam) is very similar to that 
measured simultaneously with a single dish . This implies / ~ 1. 

(4) 3C286 (Zab s ~ 0.692): IWilkinson etall 1 19791) used 609 MHz 
VLBI observations to find that 17.5 Jy of the 609 MHz flux density 
(~ 20 Jy) arises in the central 55 mas. This implie s / > 0.9. 

(5) 0458-020 (zabs ~ 2.039): iBriggs et alj 1 19891) used 608 MHz 
VLBI observations to estimate that the compact quasar core con- 
tains ~ 1.15 Jy at the redshifted Hi 21cm line frequency of ~ 
467 MHz, with an additional ~ 1 .3 Jy in two extended components, 
on scales of ~ 0".2 - 0".5 and ~ 1 - 2". They found the Hi 21cm 
optical depths measured in VLBI and single-dish studies to be in 
excellent agreement and used this to argue that the entire 2" radio 
emission is likely to be covered, i.e. / ~ 1. 

There are four DLAs with searches for Hi 21cm absorption where 
the core fraction in the background quasar is very small, imply- 
ing a low covering factor (/ < 0.1). These are at z ~ 0.437 to- 
wards 3C196, z ~ 1.3911 towards 0957+ 561 A, z~ 1.4205 to wards 
1354+258 and z ~ 0.656 towards 3C336. lBoisse et al.l J 19981) noted 
that the optical and radio sightlines towards 3C196 are clearly dif- 
ferent, due to which T s cannot be estimated. Similarly, KC03 found 
very small core fractions in 0957+561A and 1354+258, and argued 
that covering factor unc ertainties preclude a n estimate of the spin 
temperature. Conversely. ICurran et ai]( l2007l) do quote a (low) spin 
temperature for the DLA towards 3C336. However, the quasar is 
strongly lobe-dominated with the lobes extended over ~ 190 kpc, 
and only a small f raction of the flux density arises from the core 
dBridle et alj|l994l) . As in the case of 3C196, the difference be- 
tween radio and optical sightlines implies that it is not possible to 
estimate T s in this absorber. 



4 DISCUSSION 

For Hi 21cm absorption studies of DLAs towards radio-loud 
quasars, the Hi 21cm optical depth T21, Hi column density Nm 
(cm -2 ) and spin temperature T s (K) are related by the expression 



N m = 1.823 x 10 1 



(31//) / 



%2idV , 



(1) 



where / is the covering factor of the foreground DLA and the 
Hi 21cm line is assumed to be optically thin. For a DLA, where Nm 
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Table 1. Results from VLBA low-frequency imaging of quasars behind high-z DLAs. See text for details and discussion. 



is known from the Lyman-rx profile, a measurement of the Hi 21cm 
optical depth hence yields (T s /f). Estimates of the covering fac- 
tor can then be used to infer T s for the DLA. High-z DLAs have 
long been found t o have higher 7^ values than seen in the G alaxy or 
local spirals (e.g. lWolfe & Davisll 1 9791 : ICarilli et al.ll 1996] ; KC03). 
For example, KC03 found that all seven DLAs in their high-z sam- 
ple had T s > 700 K, while more than 80 % of sightlines through th e 
Milky Way and M31 have T s < 350 K dBraun & Walterbosll 1992T) . 
This suggests that high-z DLAs have far smal ler fractions of the 
cold phase of Hi than loca l spirals (see also ICarilli et al.l 1 1 9961 : 
Kanekar & Chengalur 2001). However, the T s estimates depend on 
the estimated covering fact ors. If, as argued by |CurranetalH2005h 
and|c urran & Webbl 12006), high-z DLAs have systematically low 
covering factors, the high (T s / f) measurements could arise due to 
low / values, and not due to high spin temperatures. Note that 5/7 
high-z DLAs of the KC03 sample have [T s /f] > 2000 K; extremely 
low covering factors (/ < 0.2) would then be needed to obtain spin 
temperatures in the "low" Galactic range (T s < 350 K). 

The covering factor of a DLA can be determined from VLBI 
studies at the redshifted Hi 21cm line frequency. Such studies are 
fairly common at frequencies > 1.4 GHz, implying that it is usually 
straightforward to determine T s for low-z DLAs (z < 0.5), for which 
the Hi 21cm line is redshifted to frequencies > 1 GHz. However, 
for DLAs at z > 1 .7, the line frequency is redshifted to < 500 MHz, 
where it is technically difficult to carry out VLBI studies, as the co- 
herent integra tion times are restricted by ionospheric fluctuations 
This has meant that most T s measurements in high-z 
DLAs are based on covering factors estimated either from VLBI 
studies at much higher frequencies (typically > 1. 4 GHz) or from 
the spectral index of the background quasar ( e.g. Wolfe & Davisl 
19791: IWolfe et alj|l98ll; Ide Bruvn et aUll996l: ICarilli et alj|l996l : 
Kanekar & Chengalurl 120031 ; iKanekar et alj 120061) . Prior to this 
work, there were only three DLAs a t z > 1.7 with VLB I studies at 
< 600 MHz teriggs & Wolfell 19831 ; iBriggs et alj|l989h . only one 
of which used more than two VLBI stations. 



The full sample of quasars with DLA covering factor esti- 
mates from VLBI observations now consists of 24 systems, the 
18 absorbers listed in Table[T]and six DLAs from the literature, dis- 
cussed in the last section. We have excluded the four DLAs towards 
3C196, 3C336, 0957+561 A and 1345+258 (all at z a bs < 1-5), for 
which spin temperatures are not quoted in the l iterature due to thei r 
known low covering factors [except for 3C336; rCurran et alj2007ll . 
The sample contains 14 DLAs at z > 1.5 and 10 at z < 1.5, the 
"high-z" and "low-z" sub-samples, respectively (the low-z DLAs 
are all at z < 0.9); the spatial extent of the core emission is < 500 pc 
at the DLA redshift in all cases except the z ~ 2.039 DLA to- 
wards 0458—020 (see below). No evidence for redshift evolution 
can be seen in Fig.|2tA], which plots the covering factor against 
redshift for the full sample. The median covering factor for the 
high-z sub-sample is / med ~ 0.685 ± 0.042, with / > 0.45 for all 
high-z DLAs. The low-z sample has / ~ 0.41 — 1, with a median 
value of / me d ~ 0.885±0.094. A Kolmogorov-Smirnov rank- 1 test 
yields a Gaussian probability of ~ 72% that the two sub-samples 
are drawn from different distributions, consistent (within ~ 1.1c 
significance) with the null hypothesis that they are drawn from the 
same distribution. Finally, the covering factors for the high-z sam- 
ple were all determined at frequencies < V2i cm ; these are hence 
lower limits, as the spectral index of core emission is usually in- 
verted or flat, while extended non-thermal emission typically has a 
steep spectrum, and is hence fractionally larger at lower frequen- 
cies. This renders the difference between the low-z and high-z sam- 
ples even less significant. 

The deconvolved sizes of the "core" components are < 70 mas 
for all high-z DLAs in Table [T] corresponding to spatial sizes of 
< 500 pc at the DLA redshift. Fig.|2tB] shows the spatial extent of 
the "core" emission plotted against DLA redshift; for the sources 
of Table Q] fitted by an elliptical gaussian model, the "size" is the 
geometric mean of the deconvolved major and minor axes listed 
in Column (9) of the table. In all but one of the 24 systems, the 
size is < 400 pc, which would be entirely covered even by small 
(> few kpc-sized) galaxies. The only exception is the z ~ 2.039 
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Figure 2. [A] Left panel: The covering factor of the 24 DLAs of the present VLBI sample, plotted against DLA redshift, z a bs- [B] Right panel: The spatial 
extent of the compact radio emission of the background quasars (at Zsbs) plotted versus z- d bs- The dashed line in [B] is at 400 pc. See text for discussion. 



DLA towards 0458-020, where [Briggs et al.l dl989p find that the 
entire ~ 2" radio emission is covered by the foreground DLA, im- 
plying an absorber size > 17 kpc. Note that the larger deconvolved 
sizes typically obtained at high redshifts are at least partly due to the 
fact that all high-z DLAs were observed at 327 MHz, where resid- 
ual phase errors in the data are likely to be an issue; as mentioned 
earlier, the deconvolved sizes should be treated as upper limits. 

19 of the 24 systems plotted in Fig. [2] (of which ten are at 
z > 1 - 5) have T s estima t es in the literature [e.g. IWolfe & Davisl 
1 19791 : ICarilli et al.lll996l ; iKanekar et alj|200fj : KC03]. It thus ap- 
pears that the high T s estimates in high-z DLAs are not the result of 
very low covering factors. The inferred high spin temperatures in 
high-z DLAs could th en arise due to (1) a preponderance of warm 
Hi in these systems dCarilli etal]ll996t KC03), or (2) systemati- 
cally lower Hi column densities on the radio sightlines than those 
measured towards t he optical QSO, d ue to small-scale (sub-kpc) 
structure in the Hi dWolfe et alj|2003l) . The spatial resolution of 
the present VLBA data is not sufficient to rule out the possibil- 
ity of differences between the optical and radio sightlines. How- 
ever, good agreement (within a factor of ~ 2) has been found be- 
tween Hi column densities measured along the same sightline from 
Lyman-a absorptio n and low-resolution Hi 21 cm emission studies, 
both in the Galaxy jDickev & Loc kman 199 ^) and in the z ~ 0.009 
DLA towards SBS 1549+593 dChengalur & Kanekarl 12002]) . This 
suggests that systematic large differences in Hi column density be- 
tween radio and optical sightlines in DLAs are unlikely. The high 
spin temperatures in high-z DLAs are thus more likely to be the 
result of a larger fraction of warm HI in these absorbers. 
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